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ABSTRACT 



Simulations of the magnetorotational instability (MRl) in a homogeneous 
shearing box have shown that the asymptotic strength of the magnetic field de- 
clines steeply with increasing resolution. Here I model the MRI driven dynamo as 
a large scale dynamo driven by the vertical magnetic helicity flux. This growth 
is balanced by large scale mixing driven by a secondary instability. The sat- 
urated magnetic energy density depends almost linearly on the vertical height 
of the typical eddies. The MRI can drive eddies with arbitrarily large vertical 
wavcnumber, so the eddy thickness is either set by diffusive effects, by the mag- 
netic tension of a large scale vertical field component, or by magnetic buoyancy 
effects. In homogeneous, zero magnetic flux, simulations only the first effect ap- 
plies and the saturated limit of the dynamo is determined by explicit or numerical 
diffusion. The exact result depends on the numerical details, but is consistent 
with previous work, including the claim that the saturated field energy scales as 
the gas pressure to the one quarter power (which we interpret as an artifact of 
numerical dissipation). The magnetic energy density in a homogeneous shear- 
ing box will tend to zero as the resolution of the simulation increases, but this 
has no consequences for the dynamo or for angular momentum transport in real 
accretion disks. The claim that the saturated state depends on the magnetic 
Prandtl number may also be an artifact of simulations in which microphysical 
transport coefficients set the MRI eddy thickness. Finally, the efficiency of the 
MRI dynamo is a function of the ratio of the Alfven velocity to the product of the 
pressure scale height and the local shear. As this approaches unity from below 
the dynamo reaches maximum efficiency. Farther from the disk midplane the 
Parker instability will dominate the local dynamics and the dynamo process. 
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Introduction 



For almost fifty years we liave known tliat magnetic fields embedded in a differentially 
rotating flows with a ne gative radial gradient for the angular y elocity are linearly unstable 



to radial perturbations (jVelikhov 



1959 



Chandrasekhar 



1960l ). This instability attracted 



only a modest level of interest among astrophysicists until the early 1990's when Balbus 
and Hawley (Balbus & Hawley 1991; Hawley & Balbus 1991; for a review see Balbus & 
Hawley 1998) pointed out that the magnetorotational instability (MRI) was uniquely suited 
to driving angular momentum transport in accretion disks. It is a robust linear instability, 
in the sense that the maximum growth rate is close to the local shear rate, regardless 
of the strength of the magnetic field and it feeds directly off the energy released by the 
outward transport of angular momentum. In the nonlinear limit, the angular momentum 
flux generated from this instability is proportional to the magnetic energy density. This is, 
in turn, determined by the dynamo activity generated from the turbulence. Determining 
how that dynamo operates, and the consequent mean value of in an accretion disk, then 
became a central issue in understanding ionized accretion disks. 

Early work on this problem tended to view the dynamo as a side effect of the 
turbulence, with the magnetic fleld and the kinetic energy building up together and at 



similar scales. This would natural 



approached the thermal pressure ( jBalbus fc Hawley 



ly tend to saturate as t 



1998 



le magnetic e nergy density 



Ogilvie 



20031). Alternative 



models, i nspired largely by phenomenologic al evidence from the study of cataclysmic 



variables (IMeyer fc Meyer-Hofmeister 



19831 ). looked for dynamo effects which would be 



sensitive to the geometry of the disk, i.e. the ratio of the disk height to radius. Such 
models included dynamos dr iven by waves generated by t idal interactions at large radii 



(IVishniac fc Diamond 



19921 ) and the incoherent dynamo (IVishniac &: Brandenburg 1119971 ) 



in which the large scale fleld is built up by random contributions to the electromotive force 
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from individual eddies. However, subsequent three dimensional numerical simulations of the 
MRI have failed to show evidence for a dependency on box geometry. More disturbingly, 
simulations tend to show evidence for satur ation when the magnetic energy density was 
less than 10~^ times the ambient pressure (IBrandenburg. Nordlund. Stein fc Torkelsson 



1995 



Stone et al. 



19961 ). Given phe nomenological argume nts in favor of a dimensionless 



viscosity close to 0.1 in dwarf novae (iBath fc Pringle 



19811 ) these results were somewhat 



disappointing. On the other hand, si mulations with vertica . 



saturate at substantially larger values (IMiller fc Stone 



200G 



structure have tended to 



Hawley fc Krolik 



200l|). This 



naturally raises the question of the role of magnetic buoyancy in the dynamo, both in terms 
of its effects on the MRI and the possibility that the Parker instability may contribute to 



the dynamo (ITout fc Pringle 



199^ 



Johansen fc Levin 



20081). 



More recently there has been an attempt to understand the saturation value of the 
MRI driven dynamo by returning to the simplest poss ible simulations, shearing boxes 



with no vertical gravity and no imp osed vertical field (ISano. Inutsuka. Turner fc Stone 



2004 



Fromang fc Papaloizou 



20071 ). These simulations show a steep dependence on 



numerical resolution, in the sense that at late times the magnetic energy density is inversely 
proportional to the number of grid elements. The obvious conclusion is that in the limit of 
very high resolution, or more physically, in the limit of very low viscosity and resistivity 
(and absent some externally imposed vertic al field), the MRI driven dynam o is a negligible 



effect (for a discussion of this possibility see 



Pessah. Chan &: Psaltis 



20071 ) 



In this paper we reexamine the nature of the MRI dynamo. We start from the 
assumption that it is a mean field dynamo effect, driven the magnetic helicity flux created 
by the MRI. In §2.1 we derive the dynamo growth rate. In §2.2 we find the dominant 
dissipative effect, which turns out to be driven by a secondary instability driven by vertical 
gradients in the large scale magnetic field and the angular momentum transport driven by 
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the MRI. In §2.3 we show how the typical MRI eddy thickness controls the strength of 
the magnetic field and discuss how this affects simulations and realistic accretion disks. In 
particular, we note that magnetic buoyancy provides a minimal eddy thickness for MRI 
turbulence and provides a mechanism for maintaining an MRI dynamo in a zero mean 
flux accretion disk. We briefly discuss the possible role of the Parker instability at large 
distances from the disk midplane. Finally, in §3 we summarize our results and briefly 
comment on ways to test this model. We comment, where appropriate, on the possible role 
of an imposed vertical fleld, but our treatment is aimed primarily at understanding how the 
disk dynamo functions when such environmental effects are negligible. 



2. Scaling Relations for the MRI Dynamo 



2.1. The Growth Rate for the MRI dynamo 



We consider simplest example of the MRI dynamo, a disk with no entrained magnetic 
fleld. In other words, the vertically averaged toroidal flux is zero and there is no imposed 
vertical fleld. In the limit where the vertical wavelength is much smaller than the vertical 
scale hei ghts of pressure, density or n iagnetic fleld, the linear dispersion relation has can be 



found in 



Vishniac &: Diamond 



(11992) • The presence of shear destroys the normal modes of 



the system, but they can be recovered as approximate solutions in certain limits. A detailed 



analysis can be found in 



Johnson 



(120071 ). For our purpos es the essential points a re: (1) for 



kz » one recovers the usual MRI dispersion relation (iBalbus fc Hawley 



199l|) with the 



characteristic frequency expressed as k ■ B/y/ATrp, (2) in that limit the growth rate does not 
depend on k^, (3) a locally unstable solution will grow for roughly k^/k^ e-foldings and (4) 
the characteristic radial wavenumber is of order Q/Va for all growth rates, where Q{r) is the 
angular velocity (and also the local shear) and is the Alfven speed. The MRI dispersion 
relation predicts a growth rate ~ \pik^A for small A;<^, a peak at some large fraction of f2 
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and a sharp cutoff at higher fc^. The net result is that we expect MRI driven turbulence 
to be characterized by radial and azimuthal wavenumbers of order VL/Va-, much larger 
vertical wavenumbers, and a coherence time, r, a few times Vt~^. Numerical simulations ar e 



generally consistent with these expectations (IBrandenburg. Nordlund. Stein fc Torkelsson 



1995 



Stone et al. 



19961 ). with eddies that are significantly stretched in the azimuthal 



direction, consistent with the idea that VLt is larger than one by some modest factor. 

Here we will adopt the hypothesis that the MRI drives a dynamo through the 

ic physics behind this kind of dynara o 



is set forth in 


Vishniac & Cho 


(; 


2001 


)(see 


Subramanian & Brandenburg 


200( 


3). Here 



Brandenburg fc Subramanian 



2005 



20061 ). Here we will briefly summarize the concept and the 



key results. 

Neglecting resistivity, the magnetic helicity, A ■ B, is a conserved quantity. Even 
when resistivity is not negligible the tendency of magnetic helicity to undergo an inverse 



cascade means that it is approximate 



laboratory experiments (iTaylor 



y conserved, a result that has major implications for 



19861 ). If we apply a high pass filter we can rewrite the 



magnetic helicity conservation equation in the form 

+ V • j;, + 2B • (v X b) 



-2r7(j-b), 



where /z. = (a • b) is the magnetic helicity contained in eddy scale structures, r] is the ohmic 
resistivity, and j/^ is the magnetic helicity flux contained in eddy scale structures. The 
right hand side describes the effects of dissipation and is usually negligible in astrophysical 
contexts. The third term on the left hand side of this equation gives the transfer of magnetic 
helicity between scales. This expression can be used to constrain the electromotive force, 
(v X b). The electromotive force contains a term proportional to the current helicity, (j • b), 
which is ~ k'^h in the Coulomb gauge. The effect of this term is to transfer magnetic helicity 
from small scales to large at a rate which is ~ (5^/6^) times the eddy turn over rate. When 
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this rate is faster than the dynamo growth rate we can neglect the dth term and take 



(v X b)|| ^ — V - j/,, 



(2) 



where the subscript || denotes the component parallel to the large scale magnetic field. For 
the MRI we have (w^) (b^) so we are always in th is limit. When \h is zero, or 



negligibly small, we get the phenomenon of "a quenching" (IGruzinov fc Diamond 



1994( 1. 



In the presence of differential rotation, j/^ is not zero. On dimensional grounds it has to 
be proportional to the square of the magnetic field times a transport coeffic ient. Neglecting 



constants of order unity, quasilinear estimates from 



Vishniac&Cho 



fl200lh give 



(3) 



If we follow convention and write the electromotive force as aB then this is equivalent to 



a ~ - — fir. 
Lb 



(4) 



where is the vertical scale for the large scale magnetic field. This implies a growth rate 
of 

■^^turb 



^fi ~ ^fir , 
Lb k Lb 



(5) 



where is the vertical scale length for the large scale magnetic field. For the MRI we have 
fir ~ 1, t>turb ~ B, and k^^B ~ fi. We get a dynamo growth rate of 

fi 



r 



kL 



(6) 



B 



For typical eddies we have k k^^ fi/VSi. However, since the magnetic helicity flux goes 
as it is possible, even likely, that it will be dominated by modes with kz close to VL/Va-, 
depending on the power spectrum of the MRI driven turbulence. 



A particular solution for this kind of dynamo was discussed in 



Vishniac fc Cho 



feooil), 



i.e. a dynamo in a periodic shearing box. This is appropriate for the MRI driven dynamo 
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simulations we are discussing here. Its applicability to real disks, i.e. non-periodic systems 
is less clear. This solution is characterized by a constant vertical magnetic helicity flux, 
oscillating between j/^ and large scale field components. In a real disk a necessary part 
of the dynamo process is the ejection of magnetic helicity from the t op and bottom of 



Blackman &: Field 



the d isk. However, in contrast to earlier suggestions (see, for example 
2OOII ). the ejected magnetic helicity can be contained in large scale structures or eddy scale 
structures and the net magnetic helicity of the system does not need to change. Since 
ejecting magnetic helicity bound up in large scale structures requires much less energy it is 
likely to be the preferred mode of ejection. 

The implication of this work is that efficient dynamo action is possible even in a 
periodic box provided that the dynamo process is dominated by the magnetic helicity fiux. 
There i s some work exploring the effect of b oundary conditions. We note in parti cular the 



work of 



Tobias! (ICattaneo fc Brummelll ) and 



Kapyla. Korpi fc Brandenburgl (120081 ). In both 



cases they found that periodic boundary conditions were significantly less favorable for 
the development of a large scale field than periodic boundary conditions. However, both 
simulations were in the limit of very weak large scale field where the magnetic helicity driven 
dynamo is unlikely to operate. In contrast, MRI simulations are always in the strong field 
limit. Finally, we note that recent work exploring the helicity driven dynamo in a periodic 
box with externally driven small scale turbulence has found large parts of parameter space 
with strong dynamo growth and field saturation with (5^) > (6^) (jShapovalov fc Vishniac 



20091 ). 
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2.2. Dissipation of the large scale field 

The usual estimate for the turbulent dissipation coefficient in arbitrary direction i is 

A - {vf)r ~ T^, (7) 

where fc^"^ is the correlation length in the i direction. In the quasilinear approximation we 
equate this wavenumber with the wavenumber of the dominant eddies driven by the local 
instability. For the MRl this implies a vertical diffusion coefficient which vanishes in the 
limit kz oo. In fact, the resulting turbulent diffusion rate, 

r.^i ~ (8) 

vanishes faster than the dynamo growth rate given in equation ([6]). The implication is that 
the vertical turbulent diffusion of the large scale magnetic field, driven by the MRl, is a 
negligible effect for the high resolution shearing box simulations. 

If turbulent diffusion driven by the MRl is ineffective in this case then we need to 
consider other, normally subdominant effects. Given a scalar pressure and a homogeneous 
environment there are none. However, even in the simplified environment provided by a 
homogeneous shearing box simulation of the MRl both these assumptions are violated. 
First, the stress tensor includes an anisotropic component of pViVj — bibjlAir)'^. If we 
average this over length and time scales typical of the MRl we are left with significant terms 
in the ff, 00, and f0 directions, all positive contributions of order pV^ (see, for example. 
Stone et al. 1996). Second, magnetic energy stored in large vertical wavelength modes 



provides a natural w ay to break the homq gi 



to the simulations of 



Johansen fc Youdin 



eneity of the environment. Finally, we can point 



(120081 ) which show the appearance of large scale. 



axisymmetric, slowly evolving structures in MRl turbulence in stratified and unstratified 
shearing box simulations. Evidently there is the potential for large scale vertical mixing. 
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To gain some understanding of how these structures arise we can look for hnear modes 
with slow growth rates. Shearing imposes the condition 

K > A;^-, (9) 

where lo is the magnitude of the complex frequency of our hypothetical large scale secondary 
instability. This suggests that we should concentrate on axisymmetric modes, or at least 
modes for which is negligible. Treating the MRI only as a homogeneous process that 
generates extra terms in the stress tensor we see that these terms depend on the square of 
the magnetic field. Assuming that the dynamo process is slow and that the magnetic field 
evolves only by the advection of the large scale field we can generate radial gradients in the 
extra components of the stress tensor by moving fluid elements vertically. If Pi^ = CijpV^ 
then our equations of motion are: 

dtV4> + ^Vr = -ikrCr<p2VASVA, (10) 
(dt + Drrkl)5VA = -V,d,VA, (11) 

V-'u = 0, (12) 

dtv, = —d,SP, (13) 
P 

and 

dtVr - 2VLV^ = -ikr{6P + Crr2VASVA). (14) 

Here 6P is the perturbation to the scalar pressure, SVa is the perturbation to the Alfven 
velocity, the magnetic field is assumed to be entirely in the (f> direction, and Drr is the 
radial turbulent diffusion coefficient. Given that wc expect uj << Q we can drop the dtv^ 
term in the first equation and combine the remaining equations to obtain 

{dt + D,^kl)d,v, ^ -^v^Vl (15) 
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We also note that 



(16) 



where r is the characteristic correlation time of the MRI turbulence. If we ignore D^r then 
we can find periodic traveling wave solutions to this equation with wavelengths less than 
the wavelength of the magnetic field. Multiplying both sides of equation f|T5|) by dzV^ we 
see that if f ^ oc Va then when dissipation is small dzVz increases monotonically. This is not 
the same as proving that we have an instability, but it indicates that taking Vz oc Va gives 
the best chance of an instability and that the traveling wave solutions are irrelevant. In 
order to gain further insight we take Va oc VoCos(k2;), and define two characteristic rates: 



7r = 



(17) 



and 



7d = k'^V^CrrT. 



The ratio of 7^ to 7^ can only be determined by appeal to numerical simula- 
tions, and T is difficult to estimate from the published simulations. However, it 
should be equivalent to measuring the turbulent magnetic Prandtl number, or more 
properly the Schmidt number. This can be me asured in simulations and the re- 
sults have been reported in several simulations (( 



Johansen fc Klahr 



2006 



2005 



Carballido. Stone &: Pringle 



Turner. Willacv. Bryden fc Yorke 



Johansen. Klahr fc Mee 



Carballido. Stone fc Pringle 



2006 



2005 



Fromang fc Papaloizou 



2OO6I) . The measurement varies from about 11 



2OO5I) to values of 0.8-0.9 (1 Johansen fc Klahr I l2005h . 



Johansen. Klahr &: Med (120061 ) pointed out that the value increases sharply with the applied 
vertical field and has a different dependence for the radial and vertical Schmidt numbers. 
For our purposes we ne ed the former at small ver tical field strengths, which brings us back 
to the values quoted in 



Johansen fc Klahr 



(120051). 
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If we take Vz{t = 0) = cos{kz) then we can expand the solution at all times as 

oo 

= a„(t) cos[(2n — l)/t2;]. (19) 

n=l 

In Fig. 1 we plot the ai{t) for 7^ = 0, 0.17^, and jr- The second value is a reasonably close 
approximati on to the results of Carba llido et al. whereas the last value is consistent with 



the work of 



Johansen fc Klahr 



(l2005l ). In all three cases ai increases monotonically, and 
almost exponentially. The case with no dissipation does not look like a pure mode, but the 
other two cases do. In all cases the higher order coefficients are positive, although smaller, 
and show similar behavior. These are not pure exponential modes, but the difference does 
not seem significant. It is interesting to note that a small amount of radial mixing actually 
enhances the growth of this instability. 

This linear model includes the assumption that the dynamo process is very slow. Since 
we are looking for a mixing process that will balance the dynamo, this is clearly unrealistic. 
Rather by taking this limit we are aiming for an estimate of the dissipation time scale. The 
actual dynamics of the magnetic field , with the dynamo and the secondary instability in 
statistical equilibrium, will require a more inclusive treatment, which is beyond the scope 
of this paper. 

What do we expect in the nonlinear limit? Taking krLs ~ 1, we should get circulation 
over vertical scales of order Lb on times of order L^f2/V|. This gives rise to an effective 
turbulent transport coefficient of roughly V|f2~^. This is just what we would have obtained 
if the MRI gave rise to approximately isotropic turbulence. The corresponding rms vertical 
velocity is just 

Kipple ~ LB'Jr ~ J- ^ , (20) 

which will not affect local measurements of the rms vertical velocity provided 

VI < (21) 
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Fig. 1. — The two leading Fourier coefficients for the ripphng modes as a function of time 
for varying levels of dissipation. 
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We will show below that this is the expected saturated limit for the magnetic field driven by 
the MRI. Since these modes are largely two dimensional, we can expect that they will show 
signs of inverse cascade, building coherent flows on the largest scales available. Furthermore, 
the inertial term in the z component of the force equation will produce a contribution to 
Vz which is proportional to dzV^. Since the growing mode has a large component of Vz 
proportional to Va this implies a nonlinear contribution to which is proportional d^\. 
When we substitute this into the RHS of equation f|T5l) we get a term which is independent 
of z. This cannot be balanced by the LHS of the same equation. Instead, we need to 
appeal to the neglected term, dtv^. In other words, we expect this instability to result in 
the appearance of annular bands of fast and slow rotation, supported by radial pressure 
gradients. These variations in orbital velocity will be comparable to the vertical and radial 
motions. 

Finally, we note that the slow motions driven by this rippling effect are not necessarily 
axisymmetric. Applying the condition in equation ([9]) we have 



Taking Lb to be something like the disk thickness, so that the RHS of this equation is 
roughly we see that these modes will be axisymmetric unless the disk radius is greater 
than Ha^g. In the homogeneous shearing box simulations these modes will have effectively 
no dependence. 

In terms of the evolution equation for the large scale fleld, we can include this effect by 
adding a term — -Dj-ippieJ ~ —V^Q^^J to the electromotive force. This is not quite right, 
because this creates an extra component of the electromotive force parallel to the large 
scale magnetic field. This, in turn, violates magnetic helicity conservation. In order to 
restrict the injection of magnetic helicity into the large scale field to the inverse cascade 
from the small scale field we need to subtract out the part parallel to B. In other words. 




(22) 
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the correction to the electromotive force is 

We would get the same result if we simply appealed to the force exerted by the large scale 
field and assumed t hat it was coupled to the fluid on an eddy turnover time (a heuristic 



model suggested by 



Brandenburg fc SubramanianI (l2000|)). However, despite the apparent 



agreement between the two approaches there are significant physical differences between 
them. In our approach the local shear Q enters through the coriolis force. In the ambipolar 
approximation to large scale field dynamics the quantity Q^'^ enters as the eddy turn over 
time. For the MRI these give roughly the same time scale, but this is not true in general. 
In addition, using the ambipolar approximation for the large scale field evolution in this 
way implies that we are including the pressure force exerted by the large scale field, despite 
the fact that the scalar pressure is almost unperturbed everywhere for kz very large. In 
our approach this "pressure" term is actually a description of the effects of a slow large 
scale instability and the "tension" term the result of the back reaction from small scale 
magnetic field structures arising from magnetic helicity conservation. It would appear that 
the resemblance between the two approaches is a coincidence, rooted in the properties of 
the MRI and magnetic helicity conservation. 

For an idealized a — Q dynamo, with ^ 0, \B^\ ^ Br, and no radial or azimuthal 
derivatives, we can rewrite the RHS of equation (!23l) as 

A^Hppie ~ [BldzB^f - B%Br4> - BrB^dzB^^]T. (24) 

At first glance this looks like the turbulent dissipation of -B^ is heavily suppressed, by 
a factor of ~ {Br/B^Y, but when Br and B^ are strongly correlated, as expected in a 
shearing system, we can approximate BrdzB^ ^ B^dzBr so that equation (IMl) becomes 



A^rippic ~ [BldzB^f - B^dzBr4']r. (25) 
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Then the dissipation rates for both components of the large scale magnetic field are 

A 
'2 



-di.l - ^ (26) 



Balancing this with the dynamo growth rate given in equation (E]), and using Vtr ~ 1, 
we predict a saturated magnetic field energy density of 

~ 47rp^fi2_ ^27) 

This is not quite the end, because the magnetic helicity flux becomes sharply less efficient 
at higher vertical wavenumbers, so it could easily be dominated by the smallest \kz\ modes. 
As a first approximation, we will assume that the unstable modes of the MRI are all equally 
excited. This is p hysically reasonable, and perha ps more important, it is consistent with 



the simulations of 



Fromang fc Papaloizou 



(120071 ) (see, in particular, their Fig. 6). In this 



case the magnetic helicity flux will be dominated by the largest vertical wavelength modes 
{kz ~ VL/Va) but these will contain only a fraction ~ kz/k^ of the total energy, where 
kc is the small scale cutoff to the vertical wavenumber. This implies an effective vertical 
wavenumber 



and a saturated state with 



^eff ~ ( ) , (28) 



E2 ~ AT:pLfk-'''^n\ (29) 

which we expect to apply in the asymptotic limit where /ccKa/^ ^ 1- 

Equation (l29l) is the main result of this paper, but it is physically incomplete. The 
maximum vertical wavelength, kc-, driven directly by the MRI is not a fundamental 
parameter. It is determined by effects other than the MRI itself. In real accretion disks 
kc may be set by the magnetic tension associated with some large scale vertical field. In 
stratified disks magnetic buoyancy will couple fluid motions at different heights and provide 
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a lower limit to the thickness of MRI eddies. Finally, viscosity and resistivity, either as real 
physical effects or as the consequence of limitations in numerical calculations, will impose a 
maximum vertical wavenumber. In the next section we will discuss how each of these effects 
determine the strength of large scale magnetic fields, first in numerical simulations, then in 
real accretion disks. 



2.3. The vertical wavenumber in homogeneous simulations and in accretion 

disks 

In homogeneous numerical simulations, where there is no vertical gravity, the only 
limit on the vertical wavenumber comes from the " microphysics" of the simulation, i.e. 
the numerical effects which produce an effective small scale dissipation. In the limit of 
infinite resolution, we expect —>■ oo and y| — > 0. In order to understand simulations 
with finite resolution we need to consider the scaling laws that govern in the presence 
of numerical diffusion. This requires appealing to the details of the specific numerical 
codes used to find the saturated state of the MRI in a homogeneous shearing backgroun d. 



Sano. Inutsuka. Turner fc Stone 



We w ill ex amine two recent attacks on th is problem: 
(I2OO4J ) and iFromang fc Papaloizou I (120071 ). We start with the latter because they give a 
much more detailed summary of their results. For their study Papaloizou and Fromang 
used the ZEUS code. They reported the effect of changing resolution on the distribution 
of numerical dissipation in Fourier space. While the plot of numerical dissipation versus 
wavenumber is a bit noisy, with a broad peak, they also included fits based on an equivalent 
physical dissipation. These fits, presented in their figures 9, 10, and 11, show a peak around 
k =~ 25, ~ 50, ~ 100 for vertical grid spacings of 64, 128, and 256 respectively. This allows 
us to avoid modeling the numerical dissipation and simply note that the critical value of 
kz scales almost inversely with grid spacing A in ZEUS. They also plot the reduced power 
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spectra for the kinetic and magnetic energies as a function of for all three resolutions 
in figure 6. This shows a relatively flat spectrum, in terms of the amplitude per vertical 
wavenumber, for all three resolutions, although the spectrum is fiat for a significant range 
in kz only at the highest resolution. These results confirm the notion that all the allowed 
modes in kz share equally in the MRI, and support the idea that we should weight the 
longest wavelength modes by their share of phase space, i.e. the factor Q/{kcVA), when 
estimating the magnetic helicity flux. On the other hand, equation fl29l) predicts from this 
that the magnetic energy density and angular momentum flux should then scale as A^/^. 
Fromang and Papaloizou quote a slightly steeper dependence, consistent with oc A which 
suggests that the magnetic helicity flux is dominated by the highest wavenumber modes. 
If true, this would be odd. It is more likely that once we consider the gradual flattening 
of the spectrum as the resolution is increased, the dependence on vertical wavenumber is 
consistent with our prediction. This issue could be settled by calculating the contribution 
to the magnetic helicity flux from different modes in a shearing box simulation. 

Sano et al. used a second-order Godunov-type scheme of their own devising. They did 
not report power spectra, or an effective vertical wavenumber cutoff. However, we can still 
gain some insight from their results. In their code the time step scales with the grid spacing 
divided by the magnetosonic velocity, which in this case is effectively the sound speed. It 
can be problematic to interpret numerical dissipation in terms of an effective transport 
coefficient, but for this code the most reasonable approximation is to use D = CqAcs in 
place of viscosity. Here Cq is an effective scaling for the Courant condition which decreases 
as the Courant condition becomes more stringent. The code should be able to handle MRI 
driven eddies as long as k^D is less than the MRI growth rate, or 




(30) 



In a simulaton with a fixed average density this implies kc oc P ^1^. Sano et al. reported a 
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scaling of V| oc P^/^, which imphes that in their simulations as well the saturated m agneti c 



energy density scales as ^ rather than kc We note that in a conference paper 



Sano 



( 120051 ) reported that the magnetic energy and stresses scaled linearly in the box size and 
as the square root of the grid cell dimension. This is expected from our estimate for kc in 
their work and the use of equation (|27ll rather than ( !29l) . As before, this result is broadly 
consistent with our theoretical expectations, but leaves unclear whether the deviation from 
equation fl29l) is because (/ccVa/^^)^ is not sufficiently large in these simulations, or because 
the magnetic helicity flux is dominated by the highest k^ modes. 

There is one other issue that should be raised regarding homogeneous shearing 
simulations, the role of the magnetic Prandtl number. When this is large, i.e. viscosity 
much greater than resistivity, then kc is determi ned by the viscosity a nd resistivity plays 



20071 ) show that the MRI 



no role in our model. We note that simulations (iFromang et al. 
driven dynamo operates at maximal efficiency in this regime. On the other hand, when 
the magnetic Prandtl number is small, or of order unity, then the RHS of equation ([1]) is 
no longer negligible, magnetic helicity is no longer a robustly conser v ed qu antity, and the 



efficiency of the dynamo will be reduced (as seen in iFromang et al. I (120071 )). Exactly how 
much it is reduced, and whether or not the dynamo remains viable requires us to consider 
the degree to which magnetic helicity will tend to cascade to higher wavenumbers instead 
of being transferred to very large scales. This is beyond the scope of this paper. We simply 
note that this effect will be large in any simulation in which resistivity is not much smaller 
than viscosity and in which kc is set by these constants. 

In ionized accretion disks, where the MRI is likely to be the dominant process 
for angular momentum transport, the resistivity and viscosity are very small and the 
microphysical limit on kz is so large as to be uninteresting. Consequently we can assume 
that we are in the asymptotic regime, and use equation fl29|) . There are two macrophysical 
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effects which give much smaller limits on the vertical wavenumber. First, if there is a large 
scale vertical field, possibly advected in from large radii, then is limited by the tension 
of the vertical field, i.e. 

(31) 



SO that 



/ R2 r \ 1/3 



In real disks there will be some vertical field as a consequence of the existence of large 
scale magnetic domains. Since V ■ B we have 

~ ^Br, (33) 

where is the radial magnetic domain size. This has frequently been used to argue that 
Bz ~ {H/R)Br in accretion disks. However, it is unreasonable to suppose that the typical 
magnetic domain in an accretion disk is as big the whole disk. From equation (!29|) we can 
reasonably expect the dynamo to favor magnetic domains as thick as the disk, but as long 
as Lfi > Lb all radial domain sizes are equally likely. Simple phase space considerations 
then suggest that Lji should typically be only slightly larger than Lb- Moreover, if a 
strong vertical field component can enhance the local dynamo action, this limit is strongly 
favored. All of which leads to the expectation that typically B^ is of order Br or roughly 
B(f,/ (LBkcs)- Combining equations (!28l) . (129!) . (!3T]) . and (!32l) we find that assuming that 
internally generated value of B^ determines kc results in a trivial equality. Any constraint 
imposed by B^ is scale free, i.e. independent of the parameters of the problem. Since we see 
no evidence for such a constraint in the numerical simulations, it is reasonable to conclude 
that it never comes into play. This is not to say that it is completely irrelevant. We will see 
that in a real disk the dynamo will operate with an efficiency which depends on the distance 
to the midplane. The vertical component of the large scale magnetic field constitutes a 
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vertically nonlocal effect. An efficient dynamo operating a density scale height or more 
away from the midplane can boost the efficiency of the dynamo at the midplane of the disk. 

Alternatively, we can consider the effects of vertical structure on the individual 
eddies. The MRI produces fluctuations in the local magnetic energy density of order unity, 
operating on time scales which are much longer than the sound travel time across an 
eddy. Consequently the gas pressure will vary by an amount 5P ~ and there will be 
fractional variations in the density of order / P. This will give rise to differential buoyant 
accelerations within an individual eddy of order [B'^ / P)zD? . Since these accelerations will 
act coherently for an eddy turn over time, ~ fi"^, this implies a minimal eddy thickness A2; 
of order [B'^/P)z. or 

where Lp is the pressure scale height. Using equation (l29l) for the magnetic field strength 
gives 

KLb ~ . (35) 

This will be enormous near the disk midplane, but will decrease rapidly as we move away 
from it. In terms of numerical simulations, this suggests that a background pressure scale 
height much larger than the size of the computational box will still be sufficient to remove 
the dependence on numerical resolution. For example, the highest resolution simulations 
of Fromang and Papaloizou had kcLs ~ 100, so a background pressure scale height ~ 4 
times the box size or less would have been sufficient to compete with numerical resolution 
in setting fc^- 

Equation (IMIl will be consistent with our assumption that A; ~ /c^ if 

KVa LpQ 

^ ~ » ^'^^ 
As Va approaches this limit from below the MRI eddies will become increasingly isotropic 
and the MRI dynamo will operate with increasing efficiency. Using this same expression for 
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kc, the growth rate for the MRI dynamo is 

Q Va f Va^ 



kesLs Lb \ LpVt 
The saturation level of the magnetic field is just 



(37) 



^{ATxp)[^y Lln\ (38) 

This last result shouldn't be taken too literally. It is based on a local balance between 
dissipation and the dynamo. The actual saturated state in a disk requires a global 
calculation, and must include the variation in the background density and temperature. 
The sense of this result is that the dynamo process will be much less efficient near the 
midplane, and will also depend on the diffusion of toroidal and vertical magnetic flux 
generated at larger scale heights. In the end, we expect that the magnetic energy will 
saturate at some modest fraction of the vertically averaged gas pressure, and that a more 
detailed solution requires a realistic consideration of the vertical structure of the disk. We 
note that an externally imposed vertical field will have an effect when Va,z ^i/(-^p^)r 
This means that it will have an impact first near the midplane. 

We have also assumed that the magne tic pressure is sraaller than the gas pressure. 



which is a necessary condition for the MRI (IKim fc Ostriker 
can be rewritten as 



2000). We note that this ratio 



YA = ^=(lA_\hL (39) 

This constraint is more important than Va ^ LpQ only near the disk midplane. Since 
the MRI driven dynamo is least efficient near the disk midplane, this suggests that this 
constraint is not important. 

There is one other point worth mentioning in this context. We have assumed that the 
MRI is the only instability operating in the disk. However, in a stratified disk we also expect 
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the a ppearance of the Parker instability. The hnear theory (jShu 



1974 



Foghzzo fc Tagger 



19951 ) gives us a typical growth rate of Va/Lp for azimuthal wavelengths of at least Lp. 
Consequently equation Q implies a radial wavelength comparable to Va/^- For Va ^ LpVt 
this lea ds to the expectation that the MRI will suppress the Parker instability t hrough radia l 



mixin g 



Vishniac fc Diamond 



Stone et al. 



19921 ). consistent with the numerical results of 
( 1l996l ) who first looked for signs of the Parker instability in a simulation of the MRI in 
a stratified disk. On the other hand, when Va ^ LpVL the Parker instability should be 
unaffected by the shearing of the disk. The implication is that the parameter Va/{^Lp) 
defines a continuum marked, at small values, by an MRI driven dynamo whose efficiency 
increases as Va LpQ, and at large values by a Parker instability driven dynamo. 



3. Conclusions 

In this paper we have proposed a model for the dynamo driven by the magnetorotational 
instability. It implies that the saturation strength of the magnetic field, and the attendant 
Maxwell stresses, scale as L^p^kc "^^^fl"^, i.e. slightly more steeply than linearly with box 
height and slightly less so with the MRI eddy height. The available simulations suggest 
that the dependence on each is linear. We see this in the scaling with magnetic domain 
size in the work of Sano et a 



Fromang fc Papaloizou 



The scaling with is seen most clearly in the work of 
(l2007l ). We tentatively suggest that the small difference between 
these simulations and our theoretical predictions can be blamed on the limited dynamic 
range of the simulations, i.e. kcVA not much greater than Q. The difference between the 
reported scalings with resolution seen by Sano et al. and by Fromang and Papaloizou, and 
the scaling with the gas pressure reported by Sano et al. can be explained as artifacts, due 
to the numerical methods used in the simulations. 

This model also makes specific predictions about numerical simulations of the 
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magnetorotational instability. First, the magnetic helicity flux carried by small scale 
structures should arise predominantly from eddies with vertical wavenumbers close to the 
typical radial wavenumbers. The strength of this flux should be inversely proportional to k^, 
the vertical wavenumber cutoff. Second, vertical turbulent diffusion should arise from slow 
motions on the scale of the box and should be insensitive to the vertical wavenumber cutoff. 
The radial scale of these motions should be cl ose to the box height. This seems at least 



Johansen fc Youdin 



fl2008h . The azimuthal 



qualitatively similar to the numerical results of 
scale will be larger by a factor of <y~^, i.e. usually these motions will be axisymmetric due 
to the effects of shearing. 

For real accretion disks this work implies that the behavior of the disk dynamo will be 
determined by whatever process sets the vertical correlation length of MRI turbulence. The 
most important generic effect will be magnetic buoyancy, but the presence of a vertical field 
may also be important, and can lead to an increased effective a^s? especially near the disk 
midplane. Microphysical processes that affect magnetic buoyancy may also be important. 
For example, in disks with a large radiation pressure, and the diffusion of photons from 
individual eddies, may increase the role of magnetic buoyancy, broadening eddies and 
increasing the saturation limit of the dynamo, although the same effect will increase the 
loss rate of magnetic flux from the disk Finally, viscous transport coefficients that are too 
small to suppress MRI by a factor of order unity may lead to an increased dynamo growth 
rate and a larger a^^. 

An important implication of this work is that homogeneous shearing box simulations 
cannot be used to explore the behavior of accretion disks. Without stratification and 
magnetic buoyancy the eddy thickness will be determined by numerical effects and the 
dynamo will be increasingly weak as the numerical resolution is increased. This same 
criticism can be applied to claims regarding the role of the magnetic Prandtl number, albeit 
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with less certainty. 

Finally, our results suggest a deep link between the MRI driven dynamo and the 
dynamo driven by the Parker instability. The scaling arguments presented here do not allow 
us to conclude which one accounts for the bulk of the magnetic field generation in accretion 
disks, but they do show that the two dynamos are complementary. The MRI driven dynamo 
operates at large values of Lp, i.e. near the disk midplane, but with increasing efficiency as 
we move away from the midplane. Eventually we come to a region where the MRI and the 
Parker instability are not cleanly separable from one another, and at larger heights, to a 
re gion where the Park e r inst ability drives the disk dynamo. In this context, the simulations 



of 



Johansen fc Levin 



(120081 ). which show a strong Parker instability driven dynamo in a 



magnetized accretion disk, are particularly interesting. 
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